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"<—>" ■ ABSTRACT 

O 
O 

CsJ ■ We analyze the wind generated by the great 20 year long super-Eddington outburst 

of ?7-Carinae. We show that using classical stellar atmospheres and winds theory, it is 
impossible to construct a consistent wind model in which a sufficiently small amount 
of mass, like the one observed, is shed. One expects the super-Eddington luminosity 
to drive a thick wind with a mass loss rate substantially higher than the observed one. 
The easiest way to resolve the inconsistency is if we alleviate the implicit notion that 
(NJ ' atmospheres are homogeneous. An inhomogeneous atmosphere, or "porous", allows 

more radiation to escape while exerting a smaller average force. Consequently, such 
{N| ■ an atmosphere yields a considerably lower mass loss rate for the same total luminosity. 

o 

O - 

O ■ should be revised, or at least reanalyzed carefully. 



Moreover, all the applications of the Eddington Luminosity as a strict luminosity limit 
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■ Subject headings: Radiative transfer — hydrodynamics — instabilities — stars: atmo- 

spheres — stars: individual (77 Carinae) 

> 

x. 

H , 1. Introduction 

77-Carinae is probably one of the most remarkable stellar object to have ever been documented. 
About 150 years ago, the star began a 20 year long giant eruption during which it radiated a 
supernova-like energy of roughly 3 x 10 49 ergs (Davidson & Humphreys 1997). Throughout the 
eruption it also shed some 1 — 2 Mq of material carrying approximately 6 x 10 48 ergs as kinetic 
energy (Davidson & Humphreys 1997), while expanding at a velocity of 650 km/ sec (Hiller & 
Allen 1992, Currie et al. 1996). 77-Carinae can therefore serve as a good laboratory for the study 
of atmospheres at extreme luminosity conditions. 

At first glance, it appears that the star shed a large amount of material. Indeed, the inferred 
mass loss rate during the great eruption of ~ 0.1 A^/yr is significantly larger than the mass 
loss rate inferred for the star today (^ 10 -3 A4o/yr, Davidson & Humphreys 1997 and references 
therein). However, considering that the luminosity during the great eruption is estimated to be 
significantly above the Eddington limit, we shall show that the star should have had a much higher 
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mass loss rate. In fact, it should have lost during the 20 year eruption more mass than its total 
mass, giving rise to an obvious discrepancy. 

A review of our current knowledge of 77 Car can be found in Davidson & Humphreys (1997). 
In section 2 we summarize how a wind solution for the star r] Car should be constructed. Since 
the luminosity is very high, the effects of convection must be taken into account. In section 3 we 
integrate the wind equations to show that no consistent solution for rj Car exists within the possible 
range of observed parameters. Section 4 is devoted to possible classical solutions to the discrepancy, 
showing that no such possibility exists. In section 5, we show that a porous atmosphere is a simple 
and viable solution to the wind discrepancy. 

2. Solving for the Wind 

Since the mass of 77-Car is estimated to be of order 100—120 M Q (Davidson Sz Humphreys 1997), 
the average luminosity in the great eruption was clearly super-Eddington (of the order of 5 times 
the Eddington limit). That is to say, the radiative force upwards, assuming the smallest possi- 
ble opacity (for ionized matter) given by Thomson scattering, was significantly larger than the 
gravitational pull downwards. Optically thin winds formally diverge at the Eddington limit (e.g., 
Kudritzki et al. 1989 and references therein). Consequently, a consistent wind solution requires an 
optically thick wind. We thus look for a wind in which the sonic point (which is the point at which 
the local speed of the outflow equals the speed of sound) is below the photosphere. Moreover, since 
the duration of the eruption is longer than the sound crossing time of the star by about a factor of 
50, a stationary wind appears to be a good approximation. 

In practically all super-sonic wind theories which describe super-sonic outflows from an object 
at rest, a consistent stationary solution is obtained only when the net driving force of the wind 
(excluding the pressure gradient) vanishes at the sonic point 1 . Thus, material experiencing a super- 
Eddington flux necessarily has to be above the sonic point. If most of the envelope carries a super 
Eddington flux, then no consistent stationary wind solution can be obtained and in fact, the object 
will evaporate on a dynamical time scale. In most systems however, this need not be the case. 
For example, in very hot systems (e.g., hot neutron stars during strong X-ray bursts, Quinn & 
Paczynski 1985), the opacity in the deep layers is reduced due the reduced Klein-Nishina opacity 
for Compton scattering at high temperatures. Thus, the sonic point in these objects is found where 
the temperature is high enough to reduce the opacity to the point where the flux corresponds to 
the local Eddington limit. 



lr The exception is line driven winds in which the force is explicitly a function of dv/dr which is actually an 
approximation to the line transfer equations. If we had written the proper radiation transfer equations for this case 
which only implicitly depend on dv/dr, we would have recovered that the sonic point coincides with the point at 
which the total force vanishes (cf Mihalas & Weibel Mihalas 1984 §107). Moreover, line driven winds are important 
only under optically thin conditions while we describe the optically thick part of the wind. 
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Another important effect, which should be taken into account, is convection. Deep inside the 
atmosphere, convection can carry a significant part (or almost all) of the energy flux, thus reducing 
the radiative pressure to a sub-Eddington value. In fact, as the radiative flux approaches the 
Eddington limit, convection generally arises and carries the lion share of the total energy flux (if it 
can) to keep the system at a sub-Eddington level (Joss et al. 1973). Although the total flux in the 
entire envelope (or almost all of it) can be equivalent to a super-Eddington flux, up to some depth 
below the photosphere, convection carries most of the flux so as to reduce the radiative flux alone 
into a sub-Eddington value. A consistent wind solution should therefore, have its sonic point at 
the location where the most efficient convection cannot carry enough flux any more. As we shall 
soon see, the problem in r/-Car is that this point is relatively deep within the atmosphere, where 
the density is so high that the expected mass loss is significantly higher than the observed one. 

To see this in a robust way we integrated numerically the wind equations starting from the 
photosphere inwards. The equations are those that describe optically thick spherically symmetric 
winds (Quinn & Paczynski 1985; Zytkow 1972; Kato & Hachisu 1994). The equations of mass 
conservation, momentum conservation and temperature gradient are 

4irr 2 pv = M = const (1) 

dv GM ldP„ x£r 
dr r z p dr 47rr z c 
dT _ S XP C r 2 

dr 16nr 2 caT 3y 3 X pr J ' V ' 

with standard notation. The parenthesized term in the last equation is a simple approximate 
interpolation that has the correct asymptotic limits for optical depths much larger and much smaller 
than unity (Quinn & Paczynski 1985). The last equation is the integrated form of the energy 
conservation equation. Unlike the aforementioned references, we specifically include advection by 
a maximally efficient convection. Thus, the integrated form of the energy conservation equation 
becomes 

■ fv 2 GM\ „ 

C r + Ml- + w-— \ +C conv = (4) 

C-r ~ I" £adv ~\~ l~-conv = ^-tot = £obs ~\~ £fcm,oo> (5) 

were A to t, M, C b s and liking are the total energy output of the star, the wind mass loss rate, 
the observed luminosity at infinity and the kinetic energy flux at infinity. On the other hand, 
C r , v , w, Cconvi C-adv are respectively, the local radiative luminosity, velocity, enthalpy, convective 
flux and advected flux (as internal and kinetic energies). The expression adopted for C conv is 
Airr 2 uv s where u is the internal energy per unit volume and v s is the speed of sound. By no means 
can convection be more efficient than this expression since highly dissipative shocks are unavoidable 
at higher speeds. It is likely that the maximally efficient convection is somewhat less efficient than 
this expression, but this will only aggravate the problem that we shall soon expose. Detailed 
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calculations of the wind were carried out. The calculations include the latest version of the OPAL 
opacities (Iglesias & Rogers 1996). It is found that the total opacity below the photosphere has 
comparable contributions from Thomson scattering and absorption processes. This implies that 
the modified Eddington limit, in which the Thomson opacity is replaced by the local total opacity, 
is somewhat lower than the classical Eddington limit. 

Since T e g ~ 9000°K 2 , the average luminosity implies a photospheric radius of 10 14 cm. Note 
that since it is a thick wind, the exact definition of the photosphere is ambiguous. Nevertheless, 
different definitions do not change the results by more than 10 — 20%. Knowing that the observed 
mass loss rate is roughly 0.1 Mq/ut (which gives the observed 2 Mq of shed material in 20 years, 
Davidson h Humphreys 1997), a specified flow speed at the photosphere can be translated to a 
required density. We can therefore integrate the wind equations inwards. If a consistent wind 
solution can be obtained for some value of the imposed velocity in the photosphere (which has to 
be between v s and Voo), then the integration inwards should reach a sonic point at which the total 
force on the gas vanishes. This will be attained if the convective and advective fluxes can carry a 
significant amount of the total flux so as to reduce the residual radiative flux to a sub-Eddington 
one. 



3. The Discrepancy 

We define the luminosity needed to be carried by convection and advection in order to bring 
about the vanishing of the total local force as C cr u. If enough energy is advected and convected 
then C r will be reduced to the local modified Eddington flux: 

AttcGM 

t-Edd,mod = (P) 

X 

with x the local opacity which can be larger than the Thomson opacity. Thus, from eq. (5), the 
critical advective+convective flux can be written as 

GMM 

£-crit — A-tot H ^Edd,mod- ( < ) 

r 

Figure 1 shows the fraction rj = (£ a dv + £conv) / £>crit a t the sonic point. A consistent solution 
can be found only if rj = 1 at the sonic point. 

Inspection of the figure clearly shows that the space of possible observed values does not con- 
tain a viable and consistent solution. This is of course irrespective of whether a solution from the 



2 This is the typical observed effective temperature for LBV's during outbursts (Humphreys & Davidson 1994). 
If the temperature is higher than this value, the inferred bolometric magnitude of rj Car during the eruption would 
be more negative, increasing the Eddington factor. If the temperature is lower than ~ 7000° K, the opacity at the 
photosphere and outwards rises abruptly (Davidson 1987), thus reducing the modified Eddington limit. In both cases, 
the discrepancy will be aggravated. 
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photosphere outward can or cannot be obtained. The discrepancy arises because a wind corre- 
sponding to the observed low mass loss rate necessarily has a sonic point that is not deep enough 
to have either convection or advection as an efficient mean of transporting energy. This can be seen 
from the optical depth at which the sonic point is obtained. In all cases, 1 ^ r < 300. However, 
convection is efficient only up to an optical depth of r ~ c/v s 3> 300 for p ra( i ~ p gas (Shaviv 2000b), 
or even deeper for larger radiation pressures (i.e., when close to the Eddington limit). 

4. Unfeasible Solutions to the Discrepancy 

Can the discrepancy be resolve with a classical assumption? Since the discrepancy is rather 
large, assuming the wind to emerge from an angular fraction / from the star does not relax the 
problem (it actually aggravates the problem because more material will be blown away from the 
higher luminosity regions). Another possibility that fails is having a higher velocity in the pho- 
tosphere than the one observed today for the shed material. This might be the case if the wind 
collides with previously ejected slow moving material. Even if such material did exist, the nec- 
essarily reduced mass loss rate inferred from the present day observed momentum aggravates the 
problem. 

The problem is not mitigated if we relax the assumption that the mass loss rate and the 
luminosity are assumed to be constant in time throughout the eruption. 

If one wishes to solve the problem using magnetic fields, then a solution can be found only if the 
magnetic energy density at the photosphere is significantly larger (by several orders of magnitude) 
then the equipartition value with the gas pressure. This of course seems unlikely. 

Another option is to have the distance estimate to r/-Car be three times smaller than 2300 pc. 
A shorter distance will remove r/-Car out of the cluster Trl6 of massive stars inside which it is 
observed and leave it instead roaming the inter galactic arm space. Considering the short lifetime 
of the star, just a few million years, this possibility appears as very unlikely. 

The problem can be solved if the mass of the star corresponds to a sub-Eddington luminosity. 
This proposed solution requires ??-Car to be at least a ~ 1000 Mq star. However, this suggestion is 
at variance with much lower estimates (see for instance Davidson & Humphreys 1997 and references 
therein). Nevertheless, having such a massive star is in fact not completely unrealistic and would 
have far reaching consequences if found to be true. 

5. A Viable Solution to the Discrepancy 

As the title suggests, there is a clear solution to the discrepancy. As the results show, the 
sonic point appears to be between the optical depths of ~ 1 and ~ 300. The exact value cannot 
be obtained since it requires the integration outward from the photosphere, which owing to the 
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relatively inaccurate effective temperature and therefore opacity, yields a wide range of results. If 
the mean radiative force between the point rj = 1 and the above found optical depth, is smaller 
than classically estimated, then a solution to the discrepancy can be found. Such a reduction in 
the mean radiative force is a natural result if the atmosphere is inhomogeneous. 

Shaviv (1998) has shown that in an inhomogeneous atmosphere, the effective opacity used to 
calculate the average force is reduced relative to the effective opacity used for the radiation transfer 
in a homogeneous medium. The effective opacity used for the average force should be a volume 
flux weighted average of the opacity per unit volume 3 Xv = XP- Namely, 

(XPF) 

Xeff = J^y (8) 

The effect is universal and arises in inhomogeneous systems that conduct heat or electricity. 
Extensive discussions exist in the literature under a different terminology (Isichenko 1992). The 
only requirement is therefore, that close to the Eddington limit the star develop inhomogeneities. 
The transformation from an homogeneous to an inhomogeneous atmosphere at luminosities close 
to but below Eddington luminosities, was recently found to take place generically even in Thomson 
scattering atmospheres (Shaviv 1999; Spiegel & Tao 1999; Shaviv 2000a). 

It was found that two different types of instabilities arise naturally when the luminosity ap- 
proaches the Eddington limit (Shaviv 2000a). One instability is of a phase transition into a station- 
ary nonlinear pattern of "fingers" that facilitate the escape of the radiation. The second type of 
instability allows the growth of a propagating wave, from which one expects a propagating nonlinear 
pattern to form. The two possibilities are summarized in figure 2. Both instabilities bring about a 
reduction of the average radiative force on the matter and a significant reduction of the mass loss 
rate since the sonic surface can sit near (or not much below) the photosphere. In both cases, the 
nonlinear pattern is necessarily expected to form in the region between the radius r conv at which 
rj = 1, or in other words, that C conv + C a d v is large enough to have C r ^ CEdd,modi and the photo- 
sphere. When the pattern is stationary, the rarefied regions have a larger than Eddington flux and 
the sonic surface in these regions is near r conv . On the other hand, if the pattern is propagating, 
the flux may be larger locally than the Eddington limit but the time average of the force on a 
mass element is less than Eddington. Since the instability does not occur above the photosphere, 
it should be homogeneous and hence super-Eddington with a super-sonic flow. 

Further analysis of the instabilities is needed to know which instability will dominate though 
it is more likely to be the phase transition since it is dynamically more important. 



3 When the flux is frequency dependant, a similar average should be taken in order to find the radiative force. 
However, one then takes a flux weighted mean over frequency space. 
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6. Summary 

To summarize, the super-Eddington luminosity emitted by r/-Car should have generated a much 
thicker wind with a sonic point placed significantly deeper than what can be directly inferred from 
the observations. A solution which lives in harmony with observations and theoretical modeling 
is a porous atmosphere, which allows more radiation to escape while exerting a smaller average 
force. It also means that the Eddington limit is not as destructive as one would a priori think it 
must be, even in a globally spherically symmetric case. Namely, all astrophysical analyses that 
employ the Eddington limit as a strict limit should be reconsidered carefully, even if they involve 
only unmagnetized Thomson scattering material. If r/-Carinae could have been super-Eddington 
for such a long duration without "evaporating" , other systems can display a similar behavior. 
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Fig. 1. — The fraction rj = (£ a dv + £-conv) / £>crit as a function of the photospheric velocity v p h- A 
consistent wind solution requires (a) that the velocity at the photosphere v p h satisfy: v s < v p h < v^, 
and (b) , rj(v = v s ) = 1. The thick line corresponds to the nominal observed and inferred values 
(M* = 100 M G , T eff = 9000°K, M = 0.1 M Q /yr, f C obs dt = 3 x 10 49 erg) while the additional 
lines depict the result when the values are changed to their reasonable limits (and even beyond). 
Clearly, no reasonable choice of parameters can result with a sonic point that is consistent with a 
wind solution (namely, we always find rj{v = v s ) <C 1). Basically, the discrepancy arises because 
the mass loss rate observed is too small to have the sonic point deep enough in the atmosphere 
where convection can be an efficient mean of energy transport. 
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Fig. 2. — The proposed atmospheric structure of ry Carinae during its great eruption. A homoge- 
neous atmosphere is unstable as a result of two generic instabilities that take place even in Thomson 
atmospheres when close to the Eddington limit (Shaviv 2000a). The effective opacity is therefore 
reduced (Shaviv 1998) and with it the average radiative force. The two panels describe the two 
types of possibilities for having a 'porous' atmosphere according to the characteristics of the in- 
stability that arises. An instability could produce a stationary pattern (first panel) if it originates 
from the phase transition instability and a moving pattern if it originates from the finite speed of 
light instability (second panel). See details in the text. 



